SbMYB3 transcription factor promotes root-specific flavone biosynthesis in Scutellaria baicalensis

Abstract Scutellaria baicalensis Georgi produces abundant root-specific flavones (RSFs), which provide various benefits to human health. We have elucidated the complete biosynthetic pathways of baicalein and wogonin. However, the transcriptional regulation of flavone biosynthesis in S. baicalensis remains unclear. We show that the SbMYB3 transcription factor functions as a transcriptional activator involved in the biosynthesis of RSFs in S. baicalensis. Yeast one-hybrid and transcriptional activation assays showed that SbMYB3 binds to the promoter of flavone synthase II-2 (SbFNSII-2) and enhances its transcription. In S. baicalensis hairy roots, RNAi of SbMYB3 reduced the accumulation of baicalin and wogonoside, and SbMYB3 knockout decreased the biosynthesis of baicalein, baicalin, wogonin, and wogonoside, whereas SbMYB3 overexpression enhanced the contents of baicalein, baicalin, wogonin, and wogonoside. Transcript profiling by qRT–PCR demonstrated that SbMYB3 activates SbFNSII-2 expression directly, thus leading to more abundant accumulation of RSFs. This study provides a potential target for metabolic engineering of RSFs.

Transcription factors are involved in plant growth, development, and specialized metabolism. Those controlling specialized metabolism usually function by activating or repressing transcription of 'structural' genes encoding enzymes of the regulated pathways. The R2R3-MYB family is one of the largest transcription factor families in plants, and plays important roles in plant growth and development, cell morphogenesis, primary metabolism, specialized metabolism, hormone signal transduction, and responses to environmental stresses [14][15][16][17]. MYB transcription factors comprise a DNA-binding domain of one to four MYB repeats of 50-53 amino acids, which each form helix, helix turn helix conformations. The R2R3-MYBs are the largest subfamily of MYB transcription factors in plants [11,[17][18][19][20]. The function of MYB transcription factors in crops has been studied extensively, but there are only a few related studies in medicinal plants [21]. Some studies have shown that MYB transcription factors from S. baicalensis are involved in the regulation of f lavonoid synthesis and stress responses in tobacco. Overexpression of SbMYB2 or SbMYB7 enhanced phenylpropanoid production and the tolerance of oxidative stress in transgenic tobacco, but decreased the accumulation of dicaffeoylspermidine and quercetin-3,7-O-diglucoside [22,23]. SbMYB8 regulated f lavonoid synthesis and improved plant tolerance to drought in transgenic tobacco, and electrophoretic mobility shift assays showed that this transcription factor could bind to the promoter sequence of SbCHS-2 [24]. However, the function of these SbMYBs has been studied only in heterologous hosts, and their regulatory functions as transcription factors are still largely unknown in S. baicalensis.
We aimed to identify transcription factors involved in regulation of RSF biosynthesis in S. baicalensis. We cloned the coding regions of transcription factors highly expressed in roots. Candidate SbMYB genes were screened using yeast single hybridization with the RSF-specific SbFNSII-2 promoter. Our study demonstrated that SbMYB3 can bind to the SbFNSII-2 promoter and enhances its activity, showing that it is a positive regulator of RSF biosynthesis in roots of S. baicalensis.

Results
Characterization of SbMYBs that are highly expressed in roots S. baicalensis accumulates abundant RSFs in its roots, which can be enhanced by jasmonic acid (JA) treatment. Previous studies indicated that JA promoted accumulations of baicalein, baicalin, and wogonin notably in cell suspensions or hairy roots of S. baicalensis [5,8]. To identify transcription factors that potentially regulate the synthesis of RSFs, we analyzed the transcriptomes of different organs of S. baicalensis [5]. Genes encoding six MYB proteins were found that were highly expressed in roots and could be induced by JA (Fig. 1A). We isolated fulllength cDNAs of the genes using primer pairs based on the loci identified and six SbMYBs were obtained, which contain 498-to 936-bp open reading frames (ORFs) encoding proteins of 165-311 amino acids (Supplementary Data Tables S1-S3). Among them, SbMYB3 contains a 795-bp ORF that encodes a protein of 264 amino acids with predicted molecular weight of 30.4 kDa (Supplementary Data Table S3).
Bioinformatic analysis of the six SbMYBs was carried out to identify the subgroups of R2R3-MYB transcription factors to which they belong. First, we constructed a phylogenetic tree using protein sequences of Arabidopsis thaliana MYBs and SbMYBs. Phylogenetic analysis indicated that SbMYB1, SbMYB2, SbMYB3, SbMYB4, and SbMYB6 belong to subgroup 14, 4, 20, 20, and 22 [18], respectively. SbMYB5 is most similar to AtMYB82 [25]. SbMYB3 is most similar to AtMYB62 [26], and shares 50.2% and 44.8% identity with AtMYB62 and AtMYB116 at the amino acid level, respectively (Fig. 1B). Full-length protein sequence alignment of these root-specific SbMYBs and several AtMYBs from different subgroups showed that all these SbMYBs have intact R2 and R3 domains (Supplementary Data Fig. S2) and belong to the R2R3-MYB family.
The biosynthetic pathways of baicalein and wogonin have been elucidated completely, and a gene encoding a specific isoform of f lavone synthase II (SbFNSII-2) plays a vital role in their synthesis [5,9,10]. To test whether the MYB transcription factors could bind to the SbFNSII-2 promoter, we also cloned the promoter region of −1985 to +1 upstream of the SbFNSII-2 coding region. Then, we analyzed the SbFNSII-2 promoter using the PlantCARE database [27], and found that it contains MYB-binding motifs: MYB-recognition element (MRE, AACCTAA) [28], MYB motif (CAACCA), and MYB binding cis-element (MBS, CAACTG) [29]. Moreover, other MYB-binding motifs (CACCCACCG, CACCAAA, and CACCAAA) [30] were also found in the SbFNSII-2 promoter (Supplementary Data Fig. S3). The results imply that the SbFNSII-2 promoter is a potential target of an R2R3-MYB transcription factor. We also found a Class II transposable element lying in the promoter of SbFNSII-2 (−960 to −400 upstream of the SbFNSII-2 coding region) using the Plant Repeat Database [31], and this type of transposon is a member of the hAT family. Moreover, the MYB motif lies in this transposable element (Supplementary Data Fig. S3), suggesting that this transposon may confer new root-specific regulation due to the MYB-binding site it carries.

SbMYB3 transcription factor binds to the promoter of SbFNSII-2
The yeast one-hybrid system was employed to test the potential interaction between the candidate MYBs and native SbFNSII-2 promoter. The results indicated that only the SbMYB3 transcription factor could bind to the SbFNSII-2 promoter ( Fig. 2A and Supplementary Data Fig. S4). A deletion analysis of the SbFNSII-2 promoter was carried out to determine the region where SbMYB3 binds on the SbFNSII-2 promoter. To test whether the MYB-binding motifs are responsible for the interaction between SbMYB3 and SbFNSII-2 promoter, we selected four regions of the SbFNSII-2 promoter to construct yeast bait strains: the P1 region (−1985 to −1136) without MYB-binding motifs; the P2 region (−1136 to 0) with MRE, MYB, and MBS motifs; the P3 region (−840 to 0) with MYB and MBS motifs; and the P4 region (−200 to 0) with the MBS motif (Fig. 2B). The results showed that SbMYB3 can bind to regions of P2, P3, and P4 of the SbFNSII-2 promoter, but cannot bind to the P1 region. In addition, with the decrease in the number of MYB-binding motifs in the promoter region, the strength of interaction was weakened ( Fig. 2C and D). These results suggested that the P4 region of the SbFNSII-2 promoter is adequate for the interaction with the SbMYB3 transcription factor, but the number of MYB-binding motifs determines the strength of this interaction.
In S. baicalensis, two genes, SbFNSII-1.1 and SbFNSII-1.2, encode identical proteins (SbFNSII-1) that are responsible for the synthesis of apigenin in the aerial parts of the plant [5]. The SbFNSII-2 gene lies adjacent to SbFNSII-1.2 on the same chromosome in the S. baicalensis genome, in a tail-to-tail inverted duplication, suggesting that SbFNSII-2 was produced by duplication of SbFNSII-1.2 [9]. The presence of the Class II DNA transposon in the promoter of SbFNSII-2 suggests that the duplication may have arisen as a result of aberrant transposition of a hAT transposon, Tam3, as described at the nivea locus in Antirrhinum majus [32]. We checked whether SbMYB3 could bind to the SbFNSII-1.2 promoter from −1995 to +1 upstream of the ATG start codon of SbFNSII-1. that SbMYB3 may be a transcription factor that evolved to regulate the RSF pathway. No Class II transposable element was found in the promoter region of SbFNSII-1.2.
To verify whether SbMYB3 is a transcriptional activator in plants, β-galactosidase (GUS) histochemical assays of Nicotiana tabacum hairy roots were performed. We designed three vectors: the first vector expressed only green f luorescent protein (GFP) driven by the UBI5 promoter (GFP group); the second vector expressed GFP driven by the UBI5 promoter and GUS driven by the SbFNSII-2 promoter (CK group); and the last vector expressed the same DNA regions as the second vector plus SbMYB3 driven by the CaMV35S promoter (GM group). The first and second vectors were chosen as the control groups (Fig. 3A). Then, hairy roots were induced from tobacco leaf explants transformed with Agrobacterium rhizogenes Ar1193 containing the vectors. After treatment with GUS staining solution, N. tabacum hairy roots of the GM group were stained dark blue, whereas hairy roots of the GFP group and CK groups showed no color change ( Fig. 3B and C). This work showed that SbMYB3 could trans-activate the SbFNSII-2 promoter, thus driving GUS expression. The expression patterns of SbMYB3 and GUS were consistent with the coloration results ( Fig. 3B-D). The expression of GUS was dramatically upregulated in the GM lines, and was 7-to 12-fold higher than that in the control lines (Fig. 3D). These data suggested that the SbMYB3 transcription factor can activate the SbFNSII-2 promoter, supporting its regulatory role in RSF biosynthesis.

SbMYB3 is localized in the nucleus
To determine the subcellular localization of SbMYB3, the coding region of SbMYB3 was fused in-frame to the N terminus of GFP, and the resulting construct, 35Spro::SbMYB3-GFP along the with 35Spro::GFP empty vector (EV) were transiently expressed in tobacco leaves. The results indicated that the f luorescent signal of the SbMYB3-GFP fusion protein was predominantly observed in the nucleus, and green f luorescence emitted from the SbMYB3-GFP fusion protein matched the blue f luorescence produced by DAPI (4 ,6-diamidino-2-phenylindole) staining of nuclei, whereas the signal of 35S-GFP was detected in the nucleus and cytoplasm (Supplementary Data Fig. S7), suggesting that SbMYB3 is a nucleus-localized protein.

RNAi and knockout of SbMYB3 reduced RSF biosynthesis in S. baicalensis hairy roots
RNAi of SbMYB3 was carried out in S. baicalensis hairy roots to confirm its function. Since the RNAi vector carried a red f luorescence gene (dsRED), positive hairy root lines could be screened by observation of red f luorescence. We obtained several hairy root lines with red f luorescence and screened two hairy root lines with significantly downregulated SbMYB3 expression levels. SbMYB3 expression in the RNAi lines was 42.15% and 46.58% of that in the EV line, respectively. SbMYB3 suppression decreased the levels of baicalin and wogonoside dramatically. In line 1 and line 10, baicalin levels were reduced to 1.75% and 32.90% of the control level, and wogonoside accumulation was reduced to 2.15% and 42.17% of that in the EV line, respectively. The content of baicalein was also reduced in the RNAi1 line, and was 43.90% of that in the EV line ( Fig. 4A and Supplementary Data Fig. S8A). These data suggested that SbMYB3 is a positive regulator of RSF synthesis.
Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas (CRISPR-associated protein)-mediated knockout of SbMYB3 was also performed in S. baicalensis hairy roots to confirm the role of SbMYB3. To achieve gene editing of SbMYB3, we constructed two vectors: one vector carrying dsRED driven by the CaMV35S promoter, which served as the control group; and another vector carrying dsRED driven by the CaMV35S promoter, Streptococcus pyogenes CAS9 driven by A. thaliana UBI10 promoter, and a guide RNA that targeted SbMYB3 driven by the A. thaliana U6-26 promoter (Fig. 4B). We obtained one hairy root line with mutations in both alleles (Fig. 4C). The metabolic profiles indicated that SbMYB3 knockout notably decreased the accumulations of baicalein, baicalin, wogonin, and wogonoside, which were 18.50%, 8.18%, 24%, and 16% of those in the EV line, respectively ( Fig. 4D and Supplementary Data Fig. S8B). These expression profiles were in line with patterns of f lavone accumulation, and the SbMYB3 knockout downregulated the expression of most RSF biosynthetic genes, including SbPAL2, SbPAL4, SbCLL7, SbCHS-2, SbFNSII-2, SbF6H, SbF8H, SbPFOMT5, and SbUBGAT. Their transcript levels in the knockout line were 50%, 31.6%, 7.68%, 10.9%, 28.15%, 26.24%, 8.58%, 31.89%, and 50% of those in the EV line, respectively (Fig. 4E). These results suggested that SbMYB3 knockout repressed expressions of most biosynthetic genes, thus decreasing RSF biosynthesis.

SbMYB3 overexpression enhanced RSF biosynthesis
To further verify whether SbMYB3 functions as a positive regulator for RSF biosynthesis, SbMYB3 overexpression lines (35S::SbMYB3) of S. baicalensis hairy roots were also generated. The expression profile indicated that SbMYB3 expression was upregulated in the overexpression lines, and was 2.13-to 3.97-fold higher than that in the EV line. SbMYB3 overexpression enhanced the contents of baicalein, baicalin, wogonin, and wogonoside, which were 1.77-to 3.42-, 2.25-to 4.87-, 2.32-to 7.19-, and 2.71to 5.60-fold of those in the EV line, respectively ( Fig. 5A and Supplementary Data Fig. S9). SbFNSII-2 expression was consistent with alterations of f lavone accumulation, and its level in the overexpression lines was 4.88-to 12.61-fold higher than that in the EV line, respectively. Moreover, SbF8H expression was also upregulated (Fig. 5B). These results suggested that SbMYB3 enhanced RSF biosynthesis by direct activation of SbFNSII-2 expression.

Discussion
In medicinal plants, MYB transcription factors are involved in the regulation of various specialized metabolites, such as tanshinones [33], phenolic acids [34], and artemisinin [35]. R2R3-MYBs comprise the largest MYB subfamily in plants and play vital roles in the regulation of specialized metabolism in plants, including regulation of the metabolic pathways of shikimate, terpenoid, f lavonoid, and monolignol biosynthesis [11,36,37]. Based on the conserved amino acid sequence motifs, they are classified into 28 subgroups [17]. In Arabidopsis, R2R3-MYBs in subgroup 20 are involved in stress responses (AtMYB2, AtMYB62, and AtMYB112), gibberellic acid biosynthesis (AtMYB62), jasmonate-mediated stamen maturation (AtMYB108), and regulation of anthocyanin formation (AtMYB112) [26,30,[38][39][40]. AtMYB62 has been reported to regulate phosphate starvation responses and gibberellic acid biosynthesis [26]. In this study, a new function (regulation of 4 -deoxyf lavone biosynthesis) of SbMYB3, a close homolog of AtMYB62, was elucidated in S. baicalensis. We demonstrated that the SbMYB3 transcription factor can trans-activate the SbFNSII-2 promoter, resulting in the considerable enhancement of RSF production in S. baicalensis. SbMYB3 is identical to SbMYB12 from a previous study, showing that its expression could be upregulated by treatments with abscisic acid, methyl jasmonate, and drought [41]. These results suggested that SbMYB3 is also involved in plant hormone signaling and stress responses.
Core biosynthetic genes and transcription factors involved in specialized metabolism are specifically expressed in different tissues of medicinal plants, where bioactive ingredients accumulate abundantly. RSFs are highly accumulated in S. baicalensis roots [5]. Based on the consistency of metabolic distribution and gene expression patterns, root-specific SbMYB transcription factors are candidates for regulating RSF biosynthesis (Fig. 1A). Therefore, we isolated the coding regions of six root-specific SbMYBs (Supplementary Data Tables S1-S3). SbFNSII-2 is also root-specific and plays a crucial role in RSF biosynthesis [5]. The SbFNSII-2 promoter was used to search for upstream transcription factors by yeast one-hybrid assays. We found that SbMYB3 shows strong interaction with the SbFNSII-2 promoter (Fig. 2 and Supplementary Data Fig. S4). SbFNSII-1.2 and SbFNSII-2 are separately responsible for f lavone biosynthesis in the aerial parts and roots [5], and their promoters both have MBS, MRE, and MYB motifs, but the SbFNSII-1.2 promoter could not be bound by the SbMYB3 transcription factor (Supplementary Data Figs S3, S5, and S6). These results suggested that SbMYB3 may have evolved to regulate the RSF synthesis pathway. Interestingly, the SbFNSII-2 promoter has a hAT transposon containing the MYB motif, but no Class  Significance was determined by ANOVA; * .01 < P < .05 and * * P < .01 were considered to indicate significant and highly significant levels, respectively. synthesis mainly in the roots by activating the expressions of CHS, CHI, FLAVANONE 3-HYDROXYLASE, and FLAVONOL SYNTHASE 1 (FLS1) [42]. AtMYB21 and its homologs AtMYB24 and AtMYB57 enhance f lavonol accumulation through regulation of FLS1 expression in the stamen [43]. Our work indicated that the SbMYB3 transcription factor enhanced the activity of the SbFNSII-2 promoter, which was also confirmed by the expression pattern of SbFNSII-2 under the control of SbMYB3 (Figs 3-5).
RNAi of SbMYB3 led to substantial reductions in baicalin and wogonoside accumulation, suggesting that SbMYB3 is a transcriptional activator of RSF biosynthesis (Fig. 4A). We confirmed this using CRISPR-CAS gene editing in S. baicalensis hairy roots to generate the SbMYB3 knockout line. Analysis of one line with two knockout alleles confirmed the gene silencing results and showed that knocking out SbMYB3 decreased the accumulation of baicalein, baicalin, wogonin, and wogonoside dramatically by downregulation of transcripts of SbPAL2, SbPAL4, SbCLL7, SbCHS-2, SbFNSII-2, SbF6H, SbF8H, SbPFOMT5, and SbUBGAT, confirming the role of SbMYB3 in regulating RSF biosynthesis (Fig. 4B-E). Although SbMYB3 overexpression notably enhanced RSF biosynthesis, only the expressions of SbFNSII-2 and SbF8H were upregulated, and SbFNSII-2 expression was upregulated significantly (Fig. 5). These results suggested that SbMYB3 activates RSF biosynthesis by direct control of SbFNSII-2 transcription. In Gentiana trif lora, GtMYBP3 and GtMYBP4, which belong to P1/subgroup 7 (f lavonol-specific MYBs), enhance the promoter activity of GtFNSII and regulate early f lavonoid biosynthesis in gentian f lowers [44]. CmMYB012 suppresses f lavone biosynthesis in response to high temperatures in chrysanthemum by direct regulation of CmFNS [45]. AgMYB12 binds to the AgFNS promoter and activates AgFNS expression, thus promoting apigenin biosynthesis in celery [46]. These results suggested that FNS is an important target gene which might be regulated by R2R3-MYBs involved in f lavonoid biosynthesis. The number of RSF biosynthetic genes regulated by SbMYB3 overexpression was less than that regulated by its knockout (Figs 4E and 5B). The inconsistency may be due to the fact that the upregulation of SbMYB3 expression level was much smaller than the downregulation caused by its knockout.
Taking these results together, gain or loss of function of SbMYB3 led to increased or decreased RSF accumulations, confirming that SbMYB3 is a positive regulator for RSF biosynthesis. The yeast one-hybrid assay and transcriptional activation assay demonstrated that the SbMYB3 transcription factor interacts with the SbFNSII-2 promoter and enhances its activity. Expression profiles also indicated that SbMYB3 upregulates the expression of SbFNSII-2.

Cloning and characterization of SbMYBs and the promoters of SbFNSII-1.2 and SbFNSII-2
Based on the transcriptome data of different organs and the genome of S. baicalensis [5,9], we amplified the full-length coding regions of SbMYBs, and isolated the promoters of SbFNSII-1.2 and SbFNSII-2 using specific primers (Supplementary Data  Table S1). Gene locus IDs of the isolated SbMYBs are shown in Supplementary Data Table S2. These cDNAs and promoters were constructed into vector pDONR207 and were validated by complete sequencing (Tsingke, China).

Yeast one-hybrid assay
The Matchmaker™ Gold yeast one-hybrid system (Clontech, Japan) was applied to screen the candidate transcription factors. The SbFNSII-1.2 and SbFNSII-2 promoters were separately constructed into vector pAbAi as two baits. Recombinant plasmids of proSbFNSII-1.2-pAbAi and proSbFNSII-2-pAbAi were separately digested with BbsI restriction endonuclease. Then, linearized plasmids of proSbFNSII-1.2-pAbAi and proSbFNSII-2-pAbAi were integrated into the genome of Y1HGold yeast, forming two bait-reporter yeast strains. The open reading frames of different SbMYBs were constructed into the vector pGADT7 as preys. Different recombinant plasmids of SbMYBs-pGADT7 were separately transformed into the bait strain, while the strain carrying an empty pGADT7 plasmid served as the negative control. A feasible inhibitory concentration of aureobasidin A (AbA) was applied in SD/−Leu medium to screen for transcription factors that bind to the bait. Yeast strains carrying different preys were cultured for 3 days at 28 • C. Deletion analysis of the SbFNSII-2 promoter was also performed using the above-mentioned methods.

Transcriptional activation assay
Three vectors were designed to determine if SbMYB3 is a transcriptional activator. The first vector only expresses GFP driven by the UBI5 promoter and the second vector contains GFP driven by the UBI5 promoter and GUS driven by the native SbFNSII-2 promoter, and the two vectors served as control groups. The third vector expresses the same DNA regions as the second vector plus SbMYB3 driven by the CaMV35S promoter. The three vectors were separately transformed into A. rhizogenes Ar1193. Then, successful transformants were used for induction of tobacco hairy roots. Finally, different lines of tobacco hairy roots were analyzed by GUS histochemical assays and real-time quantitative PCR (qRT-PCR).

Protein subcellular localization
The full-length cDNA (removing termination codon) of SbMYB3 was cloned into transient expression vector pBINPLUS.GFP4 and fused with the GFP gene. Plasmids of empty pBINPLUS.GFP4 and SbMYB3-pBINPLUS.GFP4 were introduced into A. rhizogenes GV3101. Positive Agrobacterium transformants were cultured in YEB medium until the OD 600 value reached 1, and centrifuged to discard the supernatant. Agrobacterium sediments were resuspended with 10 mM MgCl 2 and the OD 600 was adjusted to 1. Then, the Agrobacterium suspensions were supplemented with acetosyringone to the concentration of 100 μM and stood at room temperature for 3 hours. Agrobacterium suspensions were applied to inject young leaves of tobacco. After the injection, the tobacco was cultured in the dark for 1 day, followed by light culture for 3 days. GFP signals were observed by confocal laser microscopy (Olympus FV10i) and nuclei were stained with DAPI (Leagene, China).

Transgenic hairy root cultures
The full-length cDNA and a 300-bp fragment (non-conserved region) of SbMYB3 were introduced into pK7WG2R and pK7GWI WG2R vectors by Gateway technology, respectively. A knockout vector of SbMYB3 mediated by CRISPR was also constructed using Goldengate technology [47]. The overexpression and RNAi plasmids were transformed into A. rhizogenes A4. Due to conf lict of resistance between the knockout vector and A. rhizogenes A4, A. rhizogenes MSU440 was used to introduce the knockout plasmid. Based on a previous protocol [5], the different positive Agrobacterium transformants were applied to infect S. baicalensis leaves to induce hairy roots. Positive lines of hairy roots were confirmed by red f luorescence inspection, as the positive overexpression, RNAi, and knockout lines expressed dsRED protein. Different lines and control lines of S. baicalensis hairy roots were cultured for 50 days and collected. The hairy roots were ground into powder; one part was freeze-dried to extract f lavones and another part was used for RNA isolation.

Flavone extraction and HPLC analysis
Lyophilized hairy root powder (1.5 mg) was ultrasonically extracted with 1.5 ml of 70% methanol for 2 hours, and the powder residue was removed by centrifugation. Flavone extracts were filtered with 0.22-μm filters and then analyzed using the Agilent 1260 Infinity II HPLC system. According to a previous method [5], a 100 × 2-mm 3 μm Luna C18 (2) column was used for separation. Flavones were detected at 280 nm. Based on the retention time of standard substances and standard curves, f lavones were confirmed and measured.

Relative expression analysis by qRT-PCR
According to the protocol of the RNAprep Pure Plant Kit (TIANGEN, China), 50 mg of hairy root powder was used to isolate total RNA. cDNA was synthesized from 500 ng of RNA using the Prime-Script™ RT Master Mix (Takara, Japan). qRT-PCR was performed using TB Green ® Premix Ex Taq™ II (Takara, Japan) and specific primers (Supplementary Data Table S1). The 2 − CT method was used to calculate transcript levels of the relevant genes.